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摘 要 : 基于 1982 一 2015 年 GIMMS NDVI3g.v1 数据 ,结合 


站 点 气象 数据 ,采用 趋势 分 析 、 变 异 系数 、Hurst 指 数 及 优 


相关 分 析 等 方法 ,探讨 了 科尔沁 沙 地 植被 覆盖 的 时 空 特征 .气候 响应 及 未 来 趋势 。 结 果 表 明 :(1) 科尔沁 沙 地 34 a 
植被 覆盖 呈 绥 慢 增加 趋势 ,每 10 a 增 速 为 0.23%。 植 被 覆盖 变化 整体 上 可 分 为 "三 升 ”"1982 一 1999 年 .2000 一 2004 
年 .2008 一 2012 年 ) 和 “三 降 "(1999% 一 2000 年 .2004 一 2007 年 .2012 一 2015 年 ) 的 趋势 ,其 中 最 大 值 出 现在 1999 年 ,最 


小 值 出 现在 2009 年 。(2) 科尔沁 沙 地 植被 覆盖 格局 呈 “ 南 北 高 ,中间 低 ” 的 分 布 特征 。 以 “ 西 拉 木 伦 河 一 新 开 河 ” 为 
界线 ,北部 地 区 植被 变化 趋势 以 退化 为 主 , 南 部 地 区 以 改善 为 主 ,(3) 科尔沁 沙 地 变异 系数 西高 东 低 ,地 域 性 差异 明 
显 。 低 波动 变化 区 域 主要 分 布 在 北部 海拔 较 高 地 区 ( 占 5.52%), 其 植被 类 型 主要 为 针 阔 混交 林 。(4) 科尔沁 沙 地 植 


上 


被 变化 的 同 向 特征 强 于 反 向 特征 ,持续 退化 和 持续 改善 区 域 分 别 占 61.48% FI 37.03% 。 降 水 是 影响 研究 区 植被 变 


化 的 主要 因素 。 
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莞 漠 化 是 中 国 北方 最 主要 的 生态 环境 与 社会 
经 济 问题 之 一 "。 科 尔 沁 沙 地 在 历史 上 曾 是 水 草 丰 
美 .植被 茂盛 的 草原 ” ,但 受气 候 变 化 和 人 类 活动 的 
影响 ,逐渐 演变 成 了 中 国 面 积 最 大 的 荒漠 沙 地 ,并 
成 为 中 国 北方 沙尘暴 的 主要 物 源 之 一 * 。 植 被 是 判 
斯 生态 系统 走向 的 一 个 重要 指示 因子 ,陆地 植被 具 
有 涵养 水 源 .保持 水 土 的 作用 ,是 生态 系统 存在 的 
基础 ,因此 研究 生态 脆弱 区 植被 的 动态 变化 规律 ， 
对 于 荒漠 化 地 区 的 植物 管 护 具有 重要 意义 ”。 

植被 指数 是 衡量 地 表 植 被 状况 .指示 生态 环境 
变化 的 一 个 重要 指标 。 在 众多 植被 指数 中 , 归 一 化 
植被 指数 (Normalized Difference Vegetation Index, 
NDVI) 是 应 用 最 广泛 的 植被 指数 ,利用 NDVI 不 仅 能 
客观 反映 地 表 植 被 覆盖 状况 ,也 能 大 范围 长 时 间 地 
监测 地 表 植 被 覆盖 的 变化 5 7" 。 植 被 覆盖 受气 候 变 
化 和 人 类 活动 共同 影响 ,利用 三 大 类 植被 指数 产品 
(AVHRR, MODIS, SPOT NDVI) 开 展 不 同 尺 度 上 的 
植被 覆盖 研究 ,已 成 为 探究 全 球 生态 系统 变化 
的 一 个 重要 手段 ”。 
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Hurst 指 数 ; 偏 相关 分 析 


而 沙 地 被 公认 为 是 响应 全 球 变化 最 敏感 的 地 
区 之 一 ,土地 沙化 及 其 生态 环境 的 日 趋 恶化 对 生态 
安全 和 社会 .经济 持续 发 展 构成 严重 威胁 呈 。 科 和 尔 
沁 沙 地 是 中 国 北方 半 干 旱 农 牧 交错 带 的 典型 区 域 ， 
也 是 沙漠 化 比较 严重 的 地 区 之 一 ,分 析 响 应 全 球 
气候 背景 下 科尔沁 沙 地 的 植被 覆盖 时 空 变 化 具有 
现实 意义 。Huang 等 "" 采 用 1998 一 2006 年 SPOT 
VEGETATION NDVI 数 据 研 究 了 科尔沁 沙 地 过 去 9 
a 的 植被 覆盖 变化 ,发 现 研究 区 植被 活动 整体 呈 增 
加 趋势 , 且 增 加 幅度 较 大 。 马 龙 等 基于 2000 一 
2012 年 MODIS NDVI 数 据 对 科尔沁 沙 地 气候 与 植被 
变化 特征 及 其 定性 与 定量 响应 关系 进行 了 分 析 , 结 
果 表 明 研 究 区 植被 长 势利 东 向 西 逐 渐变 差 ;多 年 植 
被 长 势 呈 变 好 趋势 ,降水 呈 微 弱 增 加 趋势 ,气温 变 
化 趋势 不 明显 。 王 思 如 等 5 基于 自主 开发 的 水 文 
过 程 与 植被 生长 耦合 的 生态 水 文 模型 (WaVEM ) 研 
究 了 气候 变化 下 科尔沁 沙 地 植被 变化 的 响应 规律 ， 
发 现 降 水 的 年 际 变化 是 导致 植被 变化 的 主要 因素 ; 
短期 的 干旱 会 导致 植被 锐 减 ,而 干旱 结束 后 植被 能 
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够 迅速 恢复 ;多 年 连续 干旱 是 导致 科 尔 沁 植被 显著 
退化 的 主要 因素 。 魏 丰 良 等 “基于 2000 一 2009 年 
MODIS 地 表 反 射 率 产品 数据 分 析 了 科尔沁 沙 地 植 
被 覆盖 变化 及 其 与 气候 因子 的 关系 ,结果 表明 人 研究 
区 内 植被 的 生长 状况 为 东部 好 于 西部 ,南部 好 于 北 
部 ,植被 时 空 变异 性 不 大 。 李 生 勇 等 "基于 2004 一 
2013 年 的 MODIS NDVI 数 据 ,发 现 近 10a 间 科尔沁 
地 区 植被 整体 呈 增 加 趋势 ,每 年 增 速 为 3.3%。 

然而 ,针对 科尔沁 沙 地 植被 覆盖 长 期 变化 及 其 
与 气候 响应 的 研究 仍 显 不 足 , 主要 表现 在 :(1) 以 往 
的 研究 目标 过 于 单一 ,大 多 数学 者 只 针对 草地 这 一 
单一 覆盖 类 型 开展 研究 ,容易 导致 "只 见 树木 不 见 
森林 ”的 片面 认识 ;(2) 以 往 的 研究 在 气候 响应 分 析 
时 ,更 多 将 气候 变量 与 地 表 禾 被 变量 两 两 相关 分 
析 ,然而 气候 变量 之 间 是 相互 作用 的 ,多 变量 相关 
分 析 才 能 更 好 揭示 地 一 气相 互 作 用 机 制 。 鉴 于 此 ， 
本 文选 用 最 新 的 长 时 间 序 列 1982 一 2015 年 GIMMS 
NDVI3g.v1 数据 和 温度 、 降 水 数据 , 辅 以 趋势 分 析 、 
变异 系数 、Hurst 指 数 及 偏 相关 分 析 等 方法 ,对 科 尔 
沁 沙 地 地 区 的 沙 地 、 草 地 森林 和 灌 从 等 多 种 类 型 
的 植被 覆盖 变化 进行 时 空 特征 、 趋 势 特征 、 地 一 气 
耦合 特征 .可 持续 性 及 未 来 变化 趋势 等 全 面 分 析 ， 
由 在 揭示 全 球 变化 背景 下 科尔沁 沙 地 植被 覆盖 变 
化 以 及 气候 对 干旱 半 干 旱 区 植被 生长 的 影响 规律 ， 
进而 为 科尔沁 沙 地 荒漠 化 演化 趋势 评价 及 荒漠 化 
防治 工作 提供 理论 依据 。 


1 材料 与 方法 


1.1 研究 区 概况 

科尔沁 沙 地 位 于 内 蒙古 高 原 向 东北 平原 的 过 
渡 地 带 ,地理 位 置 介 于 117°49’~123°42'E、41°41'~ 
46°05'N 之 间 , 西 起 燕山 山系 的 七 老 图 山 , 东 至 松 辽 
平原 西部 , 南 以 努 鲁 儿 虎 山 为 界 , 北 接 大 兴安 岭 山 
地 南 缘 ,总 面积 为 12.51x10: ker( 图 1)。 研 究 区 
位 于 科尔沁 沙 地 腹地 ,海拔 120~800 m, HEPI E 
低 。 该 区 属于 温带 半 干 旱 大 陆 性 季风 气候 ,年 均 气 
温 5.4~6.8 C ,年 均 降水 量 343~451 mm ,年 均 风速 
3.4~4.4 ms"'。 目 前 科尔沁 沙 地 地 区 分 布 最 广 、 面 积 
最 大 的 土壤 类 型 为 风沙 质 土壤 ,其 次 是 草 甸 土 和 票 
钙 土 。 科 尔 沁 沙 地 的 原生 植被 是 以 榆树 (Ulmus 
pumila) 和 蒙古 栎 (Quercus mongolica ) 为 主 的 稀 树 草 
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图 1 研究 区 位 置 示 意图 
Fig. 1 Location ofthe Horqin Sandy Land 


原 , 但 由 于 人 为 破坏 和 沙漠 化 的 影响 ,这些 天 然 植 
被 大 部 分 已 经 退化 为 沙 生 植被 和 温带 草原 沙漠 。 
科尔沁 沙 地 处 在 中 国 北方 农 牧 交错 地 带 , 生 态 环境 
比较 脆弱 ,频繁 的 风沙 灾害 加 之 不 合理 的 人 为 活动 
导致 科尔沁 沙 地 成 为 中 国 北方 地 区 沙漠 化 最 典型 
的 地 区 之 一 ""*。 
1.2 数据 来 源 

遥感 数据 集 CIMMS NDVI3g.v1 来 源 于 美国 国 
家 航空 航天 局 (NASA) 戈 达 德 航天 中 心 ,时 间 跨 度 
H 1981—2015 4F , 28 lit] a2 HEAR W 0.083°x0.083° ,时 
间 分 辩 率 为 15 d。 该 数据 集 的 制作 经 过 了 严格 的 辐 
射 校正 .几何 校正 .去 云 及 除 坏 线 处 理 ,确保 了 数据 
质量 "”; 相 比 其 他 NDVI 数 据 集 , 该 数据 集 具 有 时 间 
序列 长 .精度 高 .误差 小 的 优点 ,在 全 球 及 区 域 植被 
变化 研究 中 具有 很 好 的 适用 性 所 。 本 研究 采用 最 
大 合成 法 (Maximum Value Composite, MVC) 获取 
1982 一 2015 年 逐 月 NDVI 数据 ,有 效 地 去 除 云 .大气 
和 太阳 高 度 角 等 产生 的 影响 ,并 进一步 采用 最 大 
合成 法 得 到 年 最 大 NDVI 数 据 集 。 

植被 类 型 数据 为 地 球 系统 科学 数据 共享 平台 
提供 的 Globalland30(https://www.webmap.cn/commres. 
do?method=globeIndex) 1:1000000 中 国 植被 类 型 图 
(2010 年 版 )。1982 一 2015 年 全 国平 均 气 温和 降水 
站 点 等 气象 数据 来 源 于 中 国 气象 科学 数据 共享 网 
(http://data.cma.cn) ,本 人 研究 利用 ArcGIS° 对 研究 区 
及 附近 100 km 范围 内 共 22 个 站 点 数据 进行 克 里 金 
(kriging) 插 值 ,从 而 获得 1982 一 2015 年 科尔沁 沙 地 
逐年 的 年 均 气 温和 降水 分 布 的 分 辨 率 为 1 km 的 机 
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格 数据 。 

1.3 研究 方法 

1.3.1 趋势 分 析 一 元 线性 回归 趋势 分 析 法 能 模拟 
每 个 栅 格 的 变化 趋势 ,反映 不 同时 期 NDVI 变化 趋 
势 的 空间 特征 。 利 用 所 有 年 份 的 数据 进行 拟 合 , 避 
免 了 研究 结果 的 随机 性 和 偶发 性 ,其 结果 可 以 更 加 
YEE BG AY JS AEE BR ZERO PB 
公式 为 : 


n 


nx $ (ix NDVI...) - [ov 
=1 =1 


i=l 


O tind = > (1) 
nx yi - by 
式 中 : 0, 为 像 元 NDVI 回 归 方 程 的 斜率 ;n 为 监测 
年 数 ; NDVI, 为 第 i; 年 的 年 最 大 NDVI 值 ; 当 0,,,> 
0 时 ,表示 随时 间 变 化 所 在 像 元 植被 指数 上 升 ,植被 
BEES MMA; RZ APE FER. 
1.3.2 变异 系数 ”变异 系数 表示 地 理 数据 的 相对 变 
化 (波动 ) 程 度 ,本 文采 用 变异 系数 来 分 析 植 被 变化 
的 稳定 性 乓 。 计 算 公 式 为 : 


heon: a- NDVI,,,) 


Se (2) 


n-1 


max 


式 中 : Cv 表示 NDVI,,,, 的 变异 系数 ;n 为 监测 年 数 ，; 
NDVI.,,, 为 第 i 年 的 年 最 大 NDVI 值 ; NDVI, 表示 
1982—2015 年 年 最 大 NDVI 的 平均 值 ; Co 值 越 大 ， 
表明 NDVI 波 动 程度 越 大 ; Cv 值 越 小 , 则 NDVI 波 动 
程度 就 越 小 。 
1.3.3 Hurst 指数 基于 重 标 极 差 (R/S) 分 析 方 法 的 
Hurst 指数 是 一 种 定量 描述 时 间 序 列 信息 长 期 依赖 
性 的 有 效 方法 ,其 能 够 有 效 预 测 时 间 序 列 数 据 相 对 
于 过 去 的 未 来 发 展 趋势 ”。 其 计算 原理 如 下 : 

给 定 一 时 间 序 列 EOY ,t=1, 2, …, 对 于 任意 正 
整数 r> 1 ,定义 均值 序列 : 


(@),=t 20 r=1,2,-- (3) 


Xe, 7)=Y(ew-(é)) 1<e<e (4) 


极 差 : 
R(7)= max X(t, 7)— min (t, T) r=1,2,-- (5) 


标准 差 : 


1 
2 


r=1,2,.. (6) 


SO= BC - 人 
考虑 比值 ROVS) = RS , 若 存在 RISx r” , 则 说 
明 时 间 序 列 EO) 存在 Hurst 现 象 ,五 值 称 为 Hurst 指 
数 , 其 值 可 在 双 对 数 坐 标 系 tnr,ln R/S) 中 用 最 小 二 
乘法 拟 合 得 到 。 
1.3.4 偏 相 关 分 析 ” 偏 相 关 分 析 用 于 研究 两 个 特定 
变量 之 间 的 相互 作用 关系 , 即 当 两 个 变量 均 与 第 三 
个 变量 存在 相关 关系 时 ,剔除 其 影响 ,而 只 分 析 另 
外 两 个 变量 之 间 的 相关 程度 2。 在 计算 偏 相关 系 
数 前 ,要 先 计 算 两 要 素 之 间 的 相关 系数 。 相 关系 数 
计算 公式 为 : 


n 


> (NDVI; a- NDVI (0-7) 
pa (7) 


| DE a NDVI.) 0 -六 


式 中 :7 为 变量 x 与 y 的 相关 系数 ;n 为 样本 数 ; 
NDVI a o y: 分 别 为 第 i 年 年 最 大 NDVI 和 年 均 气 温 
或 降水 量 ; NDVI 为 1982 一 2015 年 年 最 大 NDVI 
平均 值 ; 7 为 相应 的 时 间 段 平均 气温 或 平均 降 
水 量 。 

偏 相关 系数 计算 公式 为 : 


| iia) (ee ae 
xy xz yz 


(1 -ra)(1 | 
式 中 :7r, ,表示 变量 z 固 定 后 变量 x 和 y 的 偏 相 关 
系数 。 


2 结果 与 分 析 


2.1 植被 覆盖 时 间 变 化 特征 

1982—2015 年 科尔沁 沙 地 的 NDVI .在 0.50~ 
0.61 间 波动 (图 2) ,总 体 呈 缓慢 增加 趋势 , 与 前 人 的 
DEFER IEA BC ,每 10 a 增 速 为 0.23% , 略 高 
于 1982 一 2012 年 中 国 植被 覆盖 度 ,平均 每 10 a 增 速 
H 0.2% ,但 远 低 于 2000 一 2012 年 科尔沁 沙 地 ND- 
VI 每 10 a 整 体 增 速 为 3.4%50 和 2004 一 2013 年 科 和 尔 
泥 区 植被 覆盖 每 10 a 整体 增 速 为 33%"" 。 植 被 覆盖 
变化 整体 上 可 分 为 “三 升 "(1982 一 1999 年 ,2000 一 
2004 年 .2008 一 2012 年 ) 和 “三 降 ”"(1999 一 2000 年 、 
2004 一 2007 年 .2012 一 2015 年 ), 其 中 最 大 值 出 现在 
1999 年 ,最 小 值 出 现在 2009 年 ,这 与 王 娟 等 的 结 
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图 2 1982—2015 年 科尔沁 沙 地 植被 覆盖 时 间 变 化 趋势 
Fig. 2 Trend of vegetation coverage in Horqin Sandy 
Land during 1982-2015 


论 “1983 一 1999 年 和 2000 一 2009 年 两 个 时 段 , 科 尔 
沁 沙 地 以 南 植被 状况 两 个 时 期 均 持续 好 转 ” 以 及 魏 
丰 良 等 的 结论 “科尔沁 沙 地 植被 覆盖 时 空 变异 性 
在 时 间 上 以 2000 年 最 大 ,2005 年 最 小 ”基本 一 致 。 
1999—2015 4F [ia] , “三 降 ” 现 象 的 产生 与 同期 出 现 的 
干旱 事件 较为 吻合 

2.2 植被 覆盖 空间 变化 特征 

2.2.1 空间 分 布 特征 ”利用 1982 一 2015 年 的 逐年 
NDVI 数据 ,计算 得 到 研究 区 34 a hY NDVI n FI 
值 范围 在 0.18 ~ 0.88 之 间 , 其 空间 分 布 整 体 表现 为 : 
“南北 高 ,中间 低 ”( 图 3a)。 参 考 已 有 人 研究 ,将 科 
尔 沁 沙 地 植被 覆盖 度 分 为 5 个 等 级 : 极 低 覆盖 度 (0~ 
0.2) 、 低 覆盖 度 (0.2~0.4) P att FE (0.4~0.6) 高 覆 
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盖 度 (0.6~0.8) 、 极 高 覆盖 度 (0.8~1) WF 区 
NDVI, 的 平均 值 为 0.56, 其 中 高 值 区 占 35.49% , E 
要 分 布 在 科尔沁 沙 地 北部 海拔 450 m 以 上 的 区 域 ， 
其 主要 植被 类 型 为 阔 叶 林 、 针 叶 林 和 灌 从 ,植被 长 
势 较 好 ,此 外 ,还 有 一 部 分 分 布 在 科尔沁 南部 海拔 
较 低 、 地 势 较 平坦 的 区 域 ,主要 被 森林 和 农田 所 履 
盖 ; 中 等 覆盖 区 域 占 57.60% , 呈 大 面积 连接 成 片 分 
布 于 研究 区 的 中 部 和 南部 大 部 分 地 区 ,这 些 区 域 主 
要 被 草地 和 农田 所 覆 羡 ; 低 值 区 占 6.91% ,集中 分 布 
在 翁 牛 特 旗 东 部 以 及 奈 曼 旗 西 北部 ,这 些 地 区 主要 
为 荒漠 草地 及 裸 地 。 

2.2.2 空间 变化 趋势 ”为 检测 科尔沁 沙 地 植被 覆盖 
的 变化 趋势 ,根据 公式 (1) 逐 像 元 计算 了 1982 一 
2015 年 植被 NDVI 的 斜率 Ope ,并 对 其 进行 显著 
性 检验 (图 4)。 结 果 表 明 :(1) 研究 区 植被 覆盖 的 变 
化 趋势 具有 明显 的 空间 差异 性 ,以 “ 西 拉 木 伦 河 一 
新 开 河 "为 界线 ,北部 地 区 主要 分 布 有 农田 .草地 和 
森林 , 其 植被 覆盖 变化 以 退化 为 主 ,可 能 与 近 几 十 
年 来 气候 变异 (气温 升 高 .夏季 降水 量 减少 ”) ,呈现 
出 暖 干 化 趋势 所有 关 ; 南 部 地 区 虽然 也 分 布 有 大 量 
的 农田 和 草地 ,但 其 植被 覆盖 变化 却 以 改善 为 主 ， 
受 人 类 活动 因素 影响 明显 ,其 中 政府 实施 的 退耕 还 
林 还 草 、 禁 牧 禁 伐 及 “三 北 ” 防 护林 工程 等 一 系列 
防 沙 治 沙 工程 ,以 及 对 作物 的 科学 种 植 踢 耕地 的 科 
学 管理 .化 肥 及 农药 的 科学 使 用 等 "多 种 田间 管理 
措施 在 这 一 过 程 中 发 挥 了 重要 作用 。(2) 极 显 著 改 
善 区 域 占 研究 区 面积 的 26.49% , 主要 分 布 在 “ 西 拉 
木 伦 河 一 新 开 河 ”以 南 的 地 区 ,其 主要 植被 类 型 为 
农田 和 草地 。(3) 极 显 著 退 化 和 显著 退化 区 域 占 研 
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图 3 1982 一 2015 年 科尔沁 沙 地 平均 NDVI 空 间 分 布 及 2010 年 (Globalland 30 ) Hh 47d TAKAI 
Fig. 3 Spatial distribution of mean NDVI in Horqin Sandy Land during 1982-2015 and its land cover types in 2010 
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图 4 1982 一 2015 年 科尔沁 沙 地 NDVI 变 化 趋势 及 显著 性 
Fig.4 Trend and significance of NDVI in the Horqin Sandy Land during 1982-2015 


FE KAY 23.45% , 主要 分 布 在 “ 西 拉 木 伦 河 一 新 开 河 ” 
以 北 的 地 区 ,其 主要 植被 类 型 为 农田 、 草 地、 和 
森林 。 

223 空间 变化 稳定 性 根据 公式 (2) 计 算得 出 
1982—2015 年 科尔沁 沙 地 植被 NDVI, 的 变异 系 
数 Cv 介 于 0.019 ~0.198 之 间 ( 图 5a) ,参照 已 有 的 
研究 ,并 结合 Cv 的 大 小 以 及 本 研究 区 的 实际 情 
况 将 其 分 为 5 个 变化 等 级 : 低 波动 变化 (0.019 ~ 
0.05) 、 较 低 波 动 变化 (0.05 ~ 0.08) .中 等 波动 变化 
(0.08 ~ 0.10) 、 较 高 波动 变化 (0.10 ~ 0.13) 和 高 波动 
变化 (0.13 ~ 0.198) 。 可 以 看 出 研究 区 整体 呈现 出 
“西高 东 低 ,高 低 波动 并 存 , 较 低 波动 居多 ,地 域 性 
差异 明显 ”的 空间 分 布 特征 (图 5b)。 空 间 格 局 表现 
为 : 较 低 波动 变化 区 域 最 多 , 占 比 为 43.69% ,连接 成 
片 分 布 于 研究 区 的 东北 部 和 东南 部 ;中 等 波动 变化 
DX dak Hh 27.98% ,围绕 着 较 高 波动 区 域 ,主要 分 布 在 
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西 拉 木 伦 河 、 老 哈 河 、 新 开 河 以 及 教 来 河 的 两 侧 ; 高 
波动 变化 区 域 最 少 , 占 比 为 3.83% ,主要 分 布 在 巴林 
右 旗 中 东部 ` 奈 曼 旗 及 库伦 旗 南 部 ,其 植被 覆盖 类 
型 主要 为 草地 和 农田 (图 5b)。 进 一 步 统计 ,各 类 变 
化 中 ,草地 农田 , 裸 地 波动 性 增强 占 比 在 增加 ,和 森 
林 波 动 性 减弱 占 比 在 减少 , 灌 从 波动 性 变化 占 比 不 
明显 。 

2.3 植被 覆盖 变化 对 气候 的 响应 

在 空间 尺度 上 ,采用 偏 相关 分 析 方 法 对 科尔沁 
沙 地 植被 生长 与 气候 变化 之 间 的 响应 进行 分 析 ( 图 
6, 图 7) ,结果 表明 , NDVI， 与 降水 .气温 的 偏 相关 
系数 r。 分别 为 -0.55 ~ 0.65 ,-0.48 ~ 0.50。 

WIX NDVI 与 降水 呈正 相关 的 区 域 占 研究 
区 总 面积 的 91.85%, 其 中 具有 显著 (P<0.05 ) 和 极 显 
著 (P<0.01) 正 相关 的 区 域 占 25.65% ,主要 分 布 在 
“ 西 拉 木 伦 河 一 新 开 河 ”以 北 的 地 区 ; NDVI n 与 降 
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Fig. 5 Stability changes of NDVI in the Horqin Sandy Land during 1982-2015 
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图 6 1982 一 2015 年 科尔沁 沙 地 NDVI 与 降水 偏 相 关系 数 (a) 及 其 显著 性 (b) 
Fig.6 Partial correlation coefficient and significance between NDVI and annual average precipitation in 
Horqin Sandy Land during 1982-2015 
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图 7 1982 一 2015 年 科尔沁 沙 地 NDVI 与 气温 偏 相 关系 数 (a) 及 其 显著 性 (b) 
Fig. 7 Partial correlation coefficient and significance between NDVI and annual average temperature in 
Horqin Sandy Land during 1982-2015 


水 呈 负 相关 的 区 域 占 8.15% , EE A AG ESSE R 
东南 部 以 及 奈 曼 旗 等 地 ,其 中 具有 显著 和 极 显 著 负 
相关 的 区 域 仅 占 0.10%( 图 6)。 

研究 区 NDVI，. 与 气温 的 相关 性 整体 呈现 出 以 
“ 西 拉 木 伦 河 一 新 开 河 为 界线 , 北 负 南 正 ” 的 特征 ， 
具有 明显 的 地 域 性 差异 (图 7)。NDVI 与 气温 呈正 
相关 的 区 域 占 56.36%, 具 有 显著 和 极 显 著 正 相关 的 
区 域 占 4.32% ,主要 分 布 在 西 辽河 与 新 开 河流 经 的 
FEE .科尔沁 区 等 地 ,该 区 域 水 资源 丰富 ,受降 水 
的 影响 小 于 气温 ; NDVI， 与 气温 呈 负 相关 的 区 域 
H 43.64% ,大 多 分 布 在 研究 区 北部 ,该 地 区 蒸发 量 
大 ,年 均 气 温 较 高 而 降水 较 少 ,从 而 加 速 了 地 表 的 

进一步 将 受 温度 影响 和 受降 水 影响 主导 的 像 


元 区 分 开 , 其 中 只 受降 水 影响 的 像 元 约 占 24.50%， 
只 受气 温 影 响 的 像 元 约 占 4.08% ,同时 受气 温和 降 
水 影响 的 像 元 仅 占 1.24% ,并 日 分 布 零散 。 
2.4 植被 覆盖 未 来 变化 趋势 

科尔沁 沙 地 植被 NDVI, 的 Hurst 指 数 平均 值 
为 0.82, 其 中 大 于 0.5 的 区 域 占 98.11%, 说 明科 尔 泥 
沙 地 植被 变化 的 同 向 特征 强 于 反 向 特征 (图 8a)。 
Hurst 指数 的 高 值 区 (>0.65 ) 主 要 分 布 在 科尔沁 沙 
地 的 南部 和 北部 ,说 明 该 区 域 植被 的 未 来 变化 趋势 
与 过 去 变化 趋势 具有 很 强 的 一 致 性 ;Hurst 指 数 的 中 
值 区 (0.35<H<0.65) 主 要 分 布 在 西 拉 木 伦 河 与 乌 尔 
吉木 伦 河 的 两 侧 部 分 地 区 ,其 植被 的 未 来 变化 趋势 
具有 弱 同 向 和 弱 反 向 特征 。 

将 NDVI， 变 化 趋势 与 Hurst HAET Se HG 
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图 8 科尔沁 沙 地 NDVI 未 来 趋势 
Fig. 8 Future trend of NDVI in the Horqin Sandy Land 


合 分 析 ,进一步 揭示 植被 覆盖 的 未 来 变化 趋势 (网 
8b)。 整 体 上 看 ,科尔沁 沙 地 以 “ 西 拉 木 伦 河 一 新 开 
河 ? 为 界 ,北部 主要 呈 持 续 退 化 趋势 ,南部 呈 持 续 改 
善 趋势 。 持 续 退 化 .持续 改善 .由 改善 到 退化 及 由 
退化 到 改善 的 区 域 分 别 占 61.48% .37.03% .0.30% 及 
1.19%。 结 合 图 3b 可知, 科尔沁 沙 地 的 南部 区 域 其 
土地 覆盖 类 型 以 农田 和 人 工种 植 业 为 主 ,而 北部 多 
为 草地 和 林地 。 


3 结论 

(1) 1982 一 2015 年 科尔沁 沙 地 植被 覆盖 变化 总 
体 上 呈 缓 慢 增加 趋势 ,每 10 a 增 速 为 0.23%, 表 明科 
尔 沁 沙 地 植被 生长 发 展 状况 较 好 。 植 被 覆盖 变化 
整体 上 可 分 为 “三 升 "(1982 一 1999 年 .2000 一 2004 
年 、2008 一 2012 年 ) 和 “三 降 ”(1999 一 2000 年 、 
2004 一 2007 年 .2012 一 2015 年 ) ,其 中 最 大 值 出 现在 
1999 年 ,最 小 值 出 现在 2009 年 。 

(2) 科尔沁 沙 地 植被 覆盖 格局 呈现 出 "南北 高 ， 
中 间 低 ”的 空间 分 布 特征 ,其 中 低 值 区 占 6.91% ,该 
区 域 主要 为 荒漠 草地 和 裸 地 ;Hurst 指 数 分 析 表 明科 
尔 沁 沙 地 植被 变化 的 同 向 特征 强 于 反 向 特征 。 以 
“ 西 拉 木 伦 河 一 新 开 河 ” 为 界 , 北 部 地 区 植被 持续 退 
化 ,南部 地 区 持续 改善 , 且 持 续 退 化 和 持续 改善 区 
域 分 别 占 61.48% 和 37.03%。 

(3) 科尔沁 沙 地 植被 变异 系数 介 于 0.019~ 
0.198 之 间 ,整体 呈现 出 “西高 东 低 ,高 低 波 动 并 存 ， 
较 低 波动 居多 ,地域 性 差异 明显 ”的 特征 。 其 中 低 波 
动 变化 区 域 主要 分 布 在 北部 海拔 较 高 地 区 (5.32% ) ， 


其 植被 类 型 多 为 针 叶 林 和 落叶 阔 叶 林 ;高 波动 变化 
区 域 最 少 (3.83%) ,该 区 域 主要 为 草地 和 农田 。 

(4) NDVI，. 与 降水 和 气温 的 平均 偏 相关 系数 
分 别 为 0.23 和 0.03, 表 明 降 水 对 科尔沁 沙 地 植被 覆 


盖 的 影响 强 于 气温 。 
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Spatiotemporal variation of vegetation in the Horqin Sandy Land and its 
response to climate change from 1982-2015 
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Abstract: Based on the latest GIMMS NDVI3g.v1 data, temperature and precipitation data from 1982 to 2015, 
combined with trend analysis, variation coefficient, Hurst index and partial correlation analysis, this paper 
analyzes the spatiotemporal variation characteristics, trend characteristics, sustainability and future change trend 
of vegetation coverage in Horqin Sandy Land. The results show that: (1) the vegetation coverage of Horqin Sandy 
Land revealed a slow increasing trend during 1982-2015 (linear tendency 0.23% per 10 year). The change of 
vegetation coverage can be divided into “three ascendings” (1982-1999, 2000-2004, 2008-2012) and “three 
descendings” (1999-2000, 2004-2007, 2012-2015) on the whole, of which the maximum occurs in 1999 and the 
minimum in 2009. (2) In terms of spatial distribution characteristic, the vegetation coverage of Horqin Sandy 
Land showed a “high in the north and south and low in the middle” pattern. The low value areas accounted for 
6.91%, which were concentrated in the east of Ongniud Banner and the Northwest of Naiman Banner, and these 
areas were mainly desert grassland and bare land. (3) Taking the “Xar Moron River- Xinkai River” as the 
boundary, the vegetation in the northern area is continuously degraded, while the southern area is continuously 
improved. The areas of significant improvement accounted for 33.05%, which were mainly distributed in the 
south of “Xar Moron River- Xinkai River”. (4) Hurst index analysis showed that the same characteristics of 
vegetation coverage change were stronger than the reverse characteristics in the Horqin Sandy Land. The areas of 
continuous degradation and continuous improvement accounted for 61.48% and 37.03%, respectively. (5) The 
average partial correlation coefficients between NDVImx and precipitation and temperature were 0.23 and 0.03, 
respectively, indicating that the influence of precipitation on vegetation coverage was stronger than temperature. 
NDVI was positively correlated with precipitation, which accounted for 91.85% of the total area. 

Keywords: Horqin Sandy Land; NDVI3g.v1; spatio-temporal variation; climate response; Hurst index; partial 


correlation analysis 


